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Abstract 

By using the density functional theory method, we systematically study the influences of the 
doping of an Al atom on the electronic structures of the Mg(OOOl) surface and dissociation be- 
haviors of H2 molecules. We find that for the Al-doped surfaces, the surface relaxation around 
the doping layer changes from expansion of a clean Mg(OOOl) surface to contraction, due to the 
redistribution of electrons. After doping, the work function is enlarged, and the electronic states 
around the Fermi energy have a major distribution around the doping layer. For the dissociation 
of H2 molecules, we find that the energy barrier is enlarged for the doped surfaces. Especially, 
when the Al atom is doped at the first layer, the energy barrier is enlarged by 0.30 eV. For dif- 
ferent doping lengths, however, the dissociation energy barrier decreases slowly to the value on a 
clean Mg(OOOl) surface when the doping layer is far away from the top surface. Our results well 
describe the electronic changes after Al-doping for the Mg(OOOl) surface, and reveal some possible 
mechanisms for improving the resistance to corrosion of the Mg(OOOl) surface by doping of Al 
atoms. 
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I. INTRODUCTION 



Since hydrogen is one of the best clean fuels in the future, looking for an ideal storage 
material has become an important task for experimental and theoretical researchers. Among 
all the candidates, magnesium (Mg) based materials are widely concerned because they are 
relatively inexpensive and can hold a high weight percentage of hydrogen li So the 
interaction between hydrogen and Mg has been extensively studied js, 4, slj^,!?, ^. The main 
disadvantage that prevents the vast applications of Mg for hydrogen storage lies in that the 
corresponding hydrogenation and dehydrogenation temperatures are large . Recently, 
it is found that the aluminum (Al) doped MgAl thin films have much lower hydrogenation 
temperatures than pure Mg films ll, Q]. However, the specific reason for the improvement as 
well as the infiuence of Al doping on the electronic properties of Mg films, remains unclear. 
So, in order to advance the searches for ideal hydrogen storage materials and applications of 
hydrogen fuels, it become quite important to make clear the infiuence of Al doping on the 
electronic structure of Mg and the the interaction between hydrogen and Mg. 

Apart from the potential usages in hydrogen storage, MgAl alloys are also important 
industrial materials used in aerospace applications [lo|. It is always believed that the Al 
doping makes Mg more resist ent to being corrupted by air. So lots of efforts have been 
applied to upgrade the manufacturing techniques to dope Al in Mg materials, and now the 
technology has become quite mature. Studying the interactions between small molecules and 
Al-doped Mg surfaces can give us critical information about the mechanisms for Al doping 
to improve the resistance of Mg to corruption, and thus are very meaningful and necessary. 
Based on this background and the requirements of hydrogen storage materials, we here 
perform first principles calculations to systematically study the infiuence of Al doping on 
the electronic properties of the Mg(OOOl) surface and the adsorption behavior of hydrogen 
molecules. 

Previous studies have revealed that the Mg(OOOl) surface has both considerable s and p 



electronic states distributing around the Fermi energy, because of the sp hybridizations 



Here we further find that after doping of an Al atom, the electronic structures around the 
Fermi energy mainly distribute around the doping layer. For the dissociation of hydrogen 
molecules on the clean Mg(OOOl) surface, it has been found that the most energetically 

nnnn 

favored site is the surface bridge site p, [Tl, |8|, , with the corresponding minimum energy 
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barrier of 0.85 eV [l2|. Through our present first-principles calculations, we find that Al- 
doping will not change the most energetically favored dissociation channel. However, the 
minimum energy barrier for dissociation of hydrogen is critically dependent on the doping 
depth of Al, and the values are all larger than that on the clean Mg(OOOl) surface. The rest 
of this paper is organized as follows: In Sec. II, we describe our fist-principles calculation 
method and the models used in this paper. In Sec. Ill, we present in detail our calculated 
results, including the comparisons of the electronic structures of the Mg(OOOl) surface before 
and after Al doping, and the dissociation of hydrogen molecules on the Al-doped Mg(OOOl) 
surface. At last, the conclusion is given in Sec. IV. 



II. CALCULATIONAL METHODS 



Our calculations are performed wit 
ab-initio simulation package (VASP) 



lin density functional theory (DFT) using the Vienna 
13]. The PW91 [l^ generalized gradient approxima- 



tion and the projector- augmented wave potential [l5|] are employed to describe the exchange- 
correlation energy and the electron-ion interaction, respectively. The cutoff energy for the 
plane wave expansion is set to 250 eV, which is large enough to make the calculational 
error of the adsorption energy below 0.01 eV. The clean and Al-doped Mg(OOOl) surface 
are modeled by a slab composing of five atomic layers and a vacuum region of 20 A. The 
2x2 supercell in which each monolayer contains four Mg atoms are adopted in the study of 
the H2 adsorption. Our test calculations have shown that the 2x2 supercell is sufficiently 
large to avoid the interaction between adjacent hydrogen molecules. Integration over the 
Brillouin zone is done using the Monkhorst-Pack scheme [16] with 11 x 11 x 1 grid points. 
A Fermi broadening [17] of 0.1 eV is chosen to smear the occupation of the bands around 
the Fermi energy (Ej) by a finite-T Fermi function and extrapolating to T = K. During 
geometry optimizations, the bottom layer of the clean and Al-doped Mg(OOOl) surface is 
fixed while other Mg and Al atoms are free to relax until the forces on them are less than 
0.01 eV/A. The calculation of the potential energy surface for molecular H2 is interpolated 
to 209 points with different bond length (rfn-n) and height (/1H2) of H2 at each surface site. 
The calculated lattice constant of bulk Mg (a, c) and the bond length of a free H2 molecule 
are 3.21 A, 5.15 A and 0.75 A, respectively, in good agreement with the experimental values 
of 3.21 A, 5.20 A Q, Q and 0.74 A [20]. 
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FIG. 1: (Color online). Atomic configurations of the Mg(OOOl) (a), MgAll (b), MgA12 (c), and 
MgA13 (d) surfaces, with the corresponding supercells depicted in dashed hues. Mg and Al atoms 
are shown in grey and dark grey balls. 

III. RESULTS AND DISCUSSION 



As shown in Fig. 1, we have employed 3 different doping depths for Al in the Mg(OOOl) 
surface, with the Al atom respectively in the first (MgAll), second (MgA12) and third mono- 
layer (MgA13). For the clean Mg(OOOl) surface, our calculations reveal that the distances 
between the three topmost atomic layers are expanded from their bulk values, in accordance 



with previous theoretical js, l2l|, l22l] and experimental reports (23l]. It is because that due 
to its special surface charge redistribution, the topmost three layers of the Mg(OOOl) sur- 
face are negatively charged and hence repel each other. In fact, the surface expansion is 
a special phenomena for the Mg(OOOl) surface, which is not seen for other surfaces of Mg 



2l[ |. And instead of expansions, the surface relaxations of many metals show contractions 
^. The relaxation calculation for the clean Mg(OOOl) surface also confirms the precise of 
our methods. 

However, after the doping of an Al atom, the surface relaxation is changed. To describe 
more clearly the changes, here we define the surface relaxation as 



Adij = {dij - do)/do, 



(1) 



where dij and do are respectively the distance between the ith and jth layer of the relaxed 
surface, and the lattice spacing along the (0001) direction of bulk Mg. Table I summarizes 



TABLE I: Surface relaxations (Adi2, Ad23, and Ad^^) and work function ($) of the Mg(OOOl), 
MgAll, MgA12, and MgA13 surfaces. 

Surface relaxation 

surface $ (eV) 

Adi2 (%) Ad23 (%) Ads4 (%) 



Mg(OOOl) 


+1.95 


+0.55 


+0.82 


3.745 


MgAll 


-3.00 


-0.11 


+0.10 


3.767 


MgA12 


-2.62 


-2.60 


+0.49 


3.757 


MgA13 


+0.28 


-2.25 


-2.12 


3.757 



the calculated relaxations for the clean Mg(OOOl), MgAll, MgA12, and MgA13 surfaces. We 
can easily see from Table I that the results critically depend on the doping length of the 
Al atom. If the Al atom is doped at the first (i.e. the topmost) layer, then the distance 
between the first and second layers is contracted by 3.00 %. And if the Al atom is doped 
at the second layer, then and 6/23 respectively decreases by 2.60 % and 2.12 %. One can 
see that for the layer containing the doped Al atom, the lattice spacings with its nearest Mg 
layers are decreased a lot from that of bulk Mg. This result is confirmed by the relaxation 
results of the MgA13 surface, where g?23 and ^34 are decreased by 2.25 % and 2.12 %. By 
analyzing the atomic charges, we find that the neighboring Mg atoms always lose electrons to 
the doped Al atom. So the charging state of the doping layer and its neighboring Mg layers 
are respectively negative and positive, and the interactions between them become attracting 
forces. Thus the surface relaxation changes from expansion of the clean Mg(OOOl) surface 
to the contraction after Al doping. 

In addition, we also calculate the work function of the doped surfaces. The calculated 
values are given in Table I. We can see that the doping of Al atoms enlarges the work 
function by 0.01^0.02 eV. It means that after doping, surface electrons become harder to 
escape into the vacuum, and thus the surface's reaction activity is lowered. 

The band structures of the Mg(OOOl) and MgAll surfaces are calculated and shown in Fig. 
2, together with their electronic density of states (DOS). By comparing Figs. 2(a) and (b), 
we can see that the energy bands around the Fermi energy (Ef) change a lot after Al doping. 
Through careful wavefunction analysis, we find that the electronic states between Ei and Ef 
in Fig. 2(b) mainly distribute around the top layer. Considering that the electronic states 
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FIG. 2: (Color online). Band structures and Density of states (DOS) of the Mg(OOOl), MgAll, 
MgA12, and MgA13 surfaces. The Fermi energies are set to zero. 

around E/ are always more important in surface reactions, the adsorption and dissociation 
of small molecules on the MgAll surface are expected to be different from that on the clean 
Mg(OOOl) surface. As we will see, the dissociation of H2 molecule on the MgAll surface does 
show different characters. 

Moreover, we find that for different doping lengths, the result that the electronic states 
around Ef distribute around the doping layer still holds. As shown clearly in Fig. 3, the 
integral of the electronic states between Ei and Ef for the MgAll, MgA12, and MgA13 
surfaces respectively have major distributions around the the first, second and third layers. 
Since electrons around the first layer are more easy to take part in surface reactivities, one 
may deduce that the influence on surface reactions is weaker for Al-doping in the second and 
third layers. And we will prove that it is true for the dissociation of H2 molecules. By using 
the Bader topological method [24] for the MgAll surface, we find that the surface electron 
density around the first atomic layer also has a larger distribution around the Al atom, and 
about 3 electrons of surface Mg atoms tends to transfer to the Al atom. 

After systematical studies on the electronic states of the Mg(OOOl), MgAll, MgA12, and 
MgA13 surfaces. We build our model to study the adsorption of H2 molecules. As shown 
in Fig. 4(a), around the doped Al atom there are four different high-symmetry sites on 
the surface, respectively, the top, bridge (bri), hep and fee hollow sites. At each site, the 
adsorbed H2 molecule has three different high-symmetry orientations, respectively along the 
X (i.e., [1120]), y (i.e., [1 100]), and z (i.e., [0001]) directions. In the following, we will use top- 
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FIG. 3: (Color online). Band charge distribution around Ei? of the clean and Al-incorporation 
Mg(OOOl) surface. The energy is set to [-1.0, 0.1] eV. 



bri-x,^,z, hcp-x,^,z and fcc-x,^,z respectively to represent the 12 high-symmetry 
channels for the approaching of H2. Through systematic PES calculations, we find that 
the dissociation energy barrier of H2 along other low-symmetry channels is always larger 
than along these high-symmetry channels, similar to that observed during the dissociation 
of oxygen molecules on the Be(OOOl) surface [25I. 

By calculating the two-dimensional (2D) potential energy surface (PES) cuts, we find 
that the minimum energy path for dissociation of on the Mg(OOOl) surface is along the 
bri-^ channel, in accordance with previous reports p, Q, S]. And the corresponding energy 
barrier (AE) is 0.85 eV. After doping of an Al atom, we find that the dissociation path 
with the lowest energy barrier is still along the bri-^ channel on the MgAll, MgA12, and 
MgA13 surfaces. The calculated energy barriers, however, are all found to be larger than 
that on a clean Mg(OOOl) surface. As a typical prototype, the PES cut for a H2 molecule 
along the hii-y channel on the MgAll surface is shown in Fig. 5. And the calculated energy 
barriers for the high-symmetry adsorption channels on the four surfaces are summarized in 
Table II, together with the bond length and height of the H2 at the corresponding transition 
states. We can see from Table II that as the doping length of the Al atom is enlarged, the 
dissociation energy barrier approaches to the value of clean Mg(OOOl) surface. So the doped 
Al atom needs to be closer to the surface to make dissociation of H2 harder. 

Moreover, we can also see from Table II that for the hii-y channel, the molecular height 
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FIG. 4: (Color online), (a) The p{2 x 2) surface cell of Mg(OOOl) (clean or Al-doped) and four 
on-surface adsorption sites, (b) The sketch map showing that the molecule (with vertical or parallel 
orientation) is initially away from the surface with a hight /ih2- Mg and Al atoms are shown in 
grey and dark grey balls. The H2 molecule is depicted by red balls. 

for the transition state is smaller on the MgAll (1.04 A) than on the clean Mg(OOOl) surface 
(1.15 A). As we have already seen, the work function of the Mg(OOOl) surface is enlarged 
after Al-doping, and surface electrons become harder to escape into the vacuum. Since that 
a H2 molecule needs two more electrons to dissociate into two hydrogen atoms, the fact that 
a H2 molecule needs to be closer to the MgAll surface to dissociate is a chemical reflection 
of the change of the surface electronic states. 

For the dissociation of H2 molecules on most metal surfaces, it has been revealed that two 
diflFerent interactions exist, which are respectively the orthogonalizations between molecular 
orbitals of H2 and electronic states of metal, and the charge transfer from metal atoms to 



0, 



26 



271]. As shown in Fig. 2(b), we flnd that the Al-doping did not 



hydrogen atoms 

change the metallic properties of the Mg(OOOl) surface. So such two different interactions 
still exist for the dissociation of H2 on the MgAll surface. 

Figures 6(a) and (b) respectively show the difference electron density for the adsorption 
system of H2/MgAll before (/ih2=2.00 A) and after (/ih2=0.92 A) the transition state along 
the bri-^ channel, namely. 



Ap = p(H2 + MgAll) - p(H2) - p(MgAll)], 



(2) 
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FIG. 5: (Color online). Contour plots of the two dimensional cuts of the potential energy surfaces 
(PESs) for dissociation of H2 on the MgAU surface, as a function of the bond lengths (dn-H) and 
the heights (ha^) along the hii-y adsorption channel. 




(a) (b) 

FIG. 6: (Color online). The difference electron density for the H2/MgAll system with the height 
of H2 to be 2.00 (a) and 0.92 A (b). Mg, Al, and H atoms are shown in grey, dark grey, and red 
balls respectively. All the values are in unit of electrons per A^. 

where p(H2 + MgAll), p(H2) and p(MgAll) are respectively the electron density of the 
adsorption system, the H2 molecule and the MgAll surface. To calculate Ap, the atomic 
positions in the last two terms in Eq. 2 have been kept at those of the first term. Through 
careful wavefunction analysis, we find that at the beginning of the adsorption process, the 
molecular orbitals of H2 orthogonalize with electronic states of the MgAll surface and are 
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TABLE II: Energy barrier, bond length and height of the H2 molecule at the transition states along 
the 12 high-symmetry adsorption channels on the Mg(OOOl), MgAU, MgA12 and MgA13 surfaces. 
Energy barriers are in units of eV, while bond lengths and heights are in units of A. 
Adsorption clean Mg(OOOl) MgAll MgA12 MgA13 



channel barrier d(TS) h{TS) barrier d(TS) h{TS) barrier d{TS) h{TS) barrier d(TS) h{TS) 



bri-x 


1.60 


1.70 


1.64 


1.42 


1.54 


1.44 


1.64 


1.65 


1.52 


1.66 


1.62 


1.60 


bri-y 


0.85 


1.12 


1.15 


1.15 


1.14 


1.04 


0.91 


1.07 


1.11 


0.90 


1.08 


1.12 


top-x 


1.62 


1.58 


1.64 


1.35 


1.29 


1.38 


1.63 


1.57 


1.55 


1.61 


1.52 


1.62 


top-y 


1.63 


1.58 


1.64 


1.35 


1.29 


1.36 


1.65 


1.57 


1.54 


1.65 


1.54 


1.62 


hcp-x 


1.25 


1.34 


1.30 


1.24 


1.23 


1.19 


1.24 


1.31 


1.20 


1.21 


1.36 


1.28 


hcp-y 


1.26 


1.36 


1.31 


1.48 


1.54 


1.28 


1.25 


1.32 


1.22 


1.26 


1.36 


1.28 


fcc-x 


1.22 


2.33 


1.30 


1.22 


2.08 


1.16 


1.34 


2.27 


1.30 


1.20 


2.28 


1.28 


fcc-y 


1.23 


2.33 


1.30 


1.51 


2.48 


1.32 


1.31 


2.39 


1.27 


1.23 


2.40 


1.28 



thus broadened. This orthogonalization is very similar to that has been observed during the 
interaction of H2 with the Mg(OOOl) 12] and Al(lll) [27| surfaces. As shown in Fig. 6(a), 
there is an electron-depletion region on top of the Al atom, which is formed because electrons 
around the Al atom are repelled by the H2 bonding electrons due to the orthogonalization. 
In comparison, on the clean Mg(OOOl) surface, the electrons of Mg are repelled from the 
surface at the beginning of the adsorption of a H2 molecule. We can also see from Fig. 6(b) 
that the H2 molecule begins to get some electrons from the Mg and Al atoms, when it gets 
more closer to the MgAll surface. 

As we have already seen, extra electrons on the MgAll surface trends to distribute around 
the Al atom. Thus during the initial stage of the adsorption of a H2 molecule, the molecular 
orbitals of H2 orthogonalize only with electronic states around the Al atom, while on the 
clean Mg surface, the molecular orbitals orthogonalize with electronic states around the two 
nearest Mg atoms. Besides, the ability to lose electrons is stronger for a Mg atom than an Al 
atom. Therefore, the H2 molecule is more difficult to get electrons from the MgAll surface. 
And so the dissociation energy barrier is larger on the MgAll surface. 
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IV. SUMMARY 



In summary, the electronic structures of the Al-doped Mg(OOOl) surfaces and the dissoci- 
ation behaviors of H2 molecules on the doped surfaces are systematically studied, using the 
first-principles calculations. Our calculational results show that after Al-doping, the surface 
electronic structure is largely changed. Due to the electrons redistribution, the surface relax- 
ation changes from expansion to contraction around the doping layer. For different doping 
lengths, the work function of the Mg(OOOl) surface are all enlarged. Through careful elec- 
tronic structure analysis, we further find that the electronic states around the Fermi energy 
mainly distribute around the doping layer for Al-doped Mg(OOOl) surface. By calculating 
the potential energy surface of a H2 molecule on the Al-doped Mg(OOOl) surfaces, we find 
that the doping of an Al atom does not change the dissociation path with the lowest energy 
barrier, but enlarges the energy barrier. And for the first layer doping, the energy barrier 
is changed by 0.30 eV. From wavefunction and charge- density analysis, we reveal the two 
interaction stages during the dissociation of the H2 molecule, i.e. the first orthogonalization 
and later charge-transfer stages. We find that because more surface electrons distribute 
around the Al atom, the orthogonalization mainly happens between the H2 molecule and 
the surface Al atom. In addition, since Al atom has a weaker ability to lose electrons, the 
dissociation energy barrier for H2 is enlarged on the MgAll surface than on a clean Mg(OOOl) 
surface. Our studies provide a detailed mechanism for explaining the improvement on the 
resistance to corrosion for the Mg(OOOl) surface by doping of Al atoms. 
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